> 

o 

(N 



CX5 






X 



Semileptonic tau decays, structure functions, 
kinematics and polarisation^] The complete 
paper, including figures, is also available via 

anonymous ftp at ttp 94-27 



ttpux2.physik.uni-karlsruhe.de (129.' 

as /ttp94-27/ttp94-27.ps, or via www at 



|hLLp://LLpux2.phvsik.uni- 



karlsruhe.de/ijreijrinLs.hLml 



J. H. KUHNt 

Institut fiir Theoretische Teilchenphysik, 

Universitdt Karlsruhe, 
Kaiserstr.l2, 76128 Karlsruhe, Germany. 

and 

E. MIRKES 

Physics Department, University of Wisconsin, 
Madison, WI 53706, USA 



Abstract 



^Presented at Third Workshop on Tau Lepton Physics, Montreux, Switzerland, 19 - 22 
Sept. 1994. Work supported by BMFT contract 056 KA 93 P. 



The most general angular distribution of two or three meson final states 
from semileptonic decays r -^ TTiriy, Kniy, ttttttv, Kiririy, Kttv, KKKv, 
rjiTTTi', ... of polarized r leptons can be characterized by 16 structure 
functions. Predictions for hadronic matrix elements, based on CVC 
and chiral Lagrangians and their relations to the structure functions 
are discussed. Most of them can be determined in currently ongoing 
high statistics experiments. Emphasis of the kinematical analysis is 
firstly put on r decays in e~^e~ experiments where the neutrino escapes 
detection and the r rest frame cannot be reconstructed. Subsequently it 
is shown, how the determination of hadron tracks in double semileptonic 
events allows to fully reconstruct the r kinematics. The implications 
for the spin analysis are indicated. 



1 Structure functions 

r-decays provide an ideal for studying strong interaction physics and resonance 
properties. Detailed information on the hadronic charged current for the decay 
into three pseudoscalar mesons can be derived in particular from the study of 
angular distributions. Consider the semileptonic r-decay 

t{1,s) ~^u{l',s') + hi + h2 + h3, (1) 

into the pseudoscalar mesons hi{qi,mi) and which is governed by the matrix 
element 

M = ^{Z::)M,J^, (2) 

with G the Fermi-coupling constant. The cosine and the sine of the Cabbibo 
angle {9c) in (Q) enter for Cabibbo allowed AS = and Cabibbo suppressed 
IAS"! = 1 decays, respectively. The leptonic (M^) and hadronic (J^) currents 
are given by 

Mf, = u{l', s')-ff,{gv - gAlb)u{l, s) (3) 

and 

J'{qi.q2.q;) = {hhh\V^iO) - A^iO)\0). (4) 

V^ and A'^ are the vector and axialvector quark currents, respectively. 
The most general ansatz for the matrix element of the quark current J'^ is 
characterized by four formfactors |l| 

jM = yf^p^ + yf^p^ + ,^M^3 + VTF4 , (5) 

with 

T/M _ „M „M r)n Q{gi-g3) 
*^i — yi ^ ys ~ V q2 , 

^2 - ?2 - ft - f^ —Q-2 , .g^i 

y^^ =e'^-^^gi„g2/3g3 7 
VT =q'i+q'2+q^ =Q'- 

The formfactors Fi and F2 (-F3) originate from the axial vector hadronic cur- 
rent (vector current) and lead to spin 1 states, whereas F4 is due to the spin 
zero part of the axial current matrix element. The formfactors Fi and F2 can 
be predicted by chiral Lagrangians, supplemented by informations about res- 
onance parameters. Parametrizations of the amplitude for the Stt final states 



can be found in |l|, 0, ^ . In this case the axial fornifactors F3 is absent due to 
the G parity of the pions. (The isospin violating transition r -^ vtiuj{-^ tttt) 
would, however, give rise to such terms. See 0]) The formfactor F4 is assumed 
to be absent as a consequence of PCAC, an assumption to be tested experi- 
mentally. The 27r and Svr decay modes offer a unique tool for the study of p, p' 
resonance parameters in diffent hadronic enviroments, competing well with 
low enery e^e~ colliders with energies in the region below 1.7 GeV. As we will 
see later, the two body (p and p') resonance parameters can be determined 
in the Svr mode by taking ratios of hadronic structure functions, whereas the 
measurement of four structure functions can be used to put constraints on 
the Oi parameters. The decay modes involving different mesons (for exam- 
ple Kt^it, KKir or t^tttt) allow for axial and vector current contributions at 
the same time. Explicit parametrizations for the form factors in these decay 
modes are presented in [§], |, 0]- The vector formfactor F3 is related to the 
Wess-Zumino anomaly [Q, whereas the axial vector form factors are again pre- 
dicted by chiral lagrangians. The latter decay modes allow also for the study 
of J^*" = '" and J^*" = 1"^^ resonances which are not directly accessible in 
electron positron annihilation. 

Let us now introduce the formalism of the hadronic structure functions. 
The differential decay rate is obtained from 

dT{T -. v^ 3h) = ^ C-St:) L.^H^" dFS (7) 

where L^i, = M^^MiyY and if^"^ = J''(J'^)^ The decays are most easily 
analyzed in the hadronic rest frame where 51 + ^2 + ^3 = Q = 0. The orientation 
of the hadronic system is characterized by three Euler angles {a, (3 and 7) as 
introduced in |1], 0. In the hadronic rest frame the product of the hadronic 
and the leptonic tensors reduces to the following sum [jl| 

L^''H,, = Y,LxWx. (8) 

X 

In this system the hadronic tensor H^'^ is decomposed into 16 (real) hadronic 
structure functions Wx corresponding to 16 density matrix elements for a 
hadronic system in a spin one [Vi, V2, V^] and spin zero state [V^'] (nine of 
them originate from a pure spin one, one from a pure spin zero and six from 
interference terms). The 16 structure functions describe the dynamics of the 



three meson decay and depend only on Q^ and the Dahtz plot variables Sj. The 
factors Lx depend on the Euler angles (which determine the orientation of the 
hadronic system), on the r polarization, on the chirality parameter 7y^ of the 
ri/VT-vertex and on the total energy of the hadrons in the laboratory frame. 
This latter dependence disappears if one considers the r decay in its restframe 
where it is replaced by the emission angle of the hadron relative to the r 
spin. Analytical expressions for the 16 coefficients Lx were first presented in 
[H]. The formalism can be applied to the case, where the r rest frame cannot 
be reconstructed because of the unknown neutrino momentum and to the 
case discussed below, where the full kinematic information is available. The 
dependence of the coefficients Lx on the r polarization allows for an improved 
measurement of the r polarization at LEP [see for example 0, |10|, |Tl| and 
references therein]. 
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Figure 1: Ratio of the spin one hadronic structure 

functions wc/wajWd/wa^weIwa (from top to bottom) for r -^ uttiitt as a 
function of Q^. 



The hadronic structure functions Wx on the other hand contain the full 
dynamics of the hadronic decay. The measurement of these structure functions 
provides a unique tool for low enery hadronic physics. They can be calculated 



from the components of the hadronic current J^ and expressed in terms of the 
form factors Fj. For brevity only the results for the pure spin one state are 
listed. 



Wa = 


{xl + xD Fi\^ + {xl + xl) F2P 




+ 2{xiX2- xl)Re {FiF;) , 


Wb = 


xl F3P 


Wc = 


[xl-xl)\F,^ + {xl-xl)\F,^ 



+ 2(xiX2 + x2)Re(FiF*) 

Wd = 2[xiX3|Fi|2-X2X3|F2|2 

+ X3{X2 - xi)Re{FiF;)] 
We = -2x3(xi + X2)Im(FiF2*) (9) 

Wf = 2x4[xilm{FiF;) + X2lm{F2F;)] 
Wg = -2x,[x,Re{F,F;)+X2Re{F2F;)]] 
Wh = 2x3X4 [Im (F1F3*) - Im (F2F3*)] 

Wj = -2x,x,[Re{F,F;)-Re{F2F;)]. 

The remaining structure functions originating from a possible (small) contri- 
bution of a spin zero state are presented in [1]. The variables Xj are defined 

by xi = V^ = <?! — Q'f, X2 = V2 = ^2 ~ ^f; ^3 — ^1 = Qi = ~Q2y ^4 = ^i = 
A/Q^xsgf , where qf (gf) denotes the x (y) component of the momentum of 
meson i in the hadronic rest frame. They can easily be expressed in terms of 
Si, S2 and S3 fl], H. The structure functions can be extracted by taking suit- 
able moments with respect to an appropriate product of sine or cosine of two 
Euler angles. An alternative method to extract the structure functions has 
been suggested in [|1^ where a direct fit to the expressions (H) was performed. 



As an example, we will now present numerical results for the non vanishing 
structure functions Wa, Wc, Wd and We in the Svr decay mode. Figure 1 shows 
predictions for the structure function ratios wc/wa, wd/wa and we/wa as a 
function of Q^, where we have integrated over the Dalitz plot variables Sj [The 
integrated structure functions are denoted by lower case letter wx] ■ The results 



are based on the same parametrization of the forinfactors as those used in ||l|]. 
Although information on the resonance parameters in the two body decays 
is lost by integrating over Si and S2, interesting structures Q are observed 
nevertheless. One observes that all normalized structure functions are sizable. 
wc/wa approaches its maximal value 1 for small Q^. The "time reversal" 
invariant ratios wcIwa and wd/wa reach up to unity, the "non-invariant" 
ratio we/wa which is sensitive towards phases is significantly smaller. Note, 
that the dependence on the ai mass and width parameters cancel in the ratio 
wx/wa in fig. 1. 
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Figure 2: Spin one hadronic structure functions wa,Wc,Wd,'We (from top to 
bottom) for r —>■ z/tttttt as a function of Q^. Results are shown for two sets of 
ai parameters: nia^ = 1.251 GeV, Tai = 0.599 GeV (solid) and niai = 1.251 



GeV, r^, = 0.550 GeV (dashed) 



On the other hand, the Q^ distributions of the structure functions wa,c,d,e 
presented in fig. 2 are very sensitive to the Oi parameters. As an example, fig. 
2 shows predictions for the structure functions, where two different values for 



-'^More constraints on the two body resonances can be obtained by analyzing the fuU 
dependence on Q^ and Si, which should be accessible with the present high statistic 
experiments. 



the ai width have been used. Therefore, the ratios in fig. 1 can be used to 
fix the model dependence in the two body resonances, whereas the structure 
functions themselves impose rigid contraints on the Oi parameters. Through 
measurements of the Wx it is therefore possible to determine the amplitudes 
in much more detail than through rate measurements alone. 

The technique of structure functions also allow for a model independent 
search for a spin zero component in the hadronic current. Such a contribution 
would lead to the six forementioned additional structure functions originating 
from the interference with the (large) spin one contributions|l[|. 

The analysis of angular distributions of the hadronic final states allows 
to determine not only the properties of the hadronic current, it is also an 
important tool for the determination of the r polarization and the structure of 
the TuW coupling. This holds true for single pion final states, but is equally 
valid for decays into two|]12|, |13|, |n| or three |T2|, |1^, |1], ^ pions, if the angular 



distributions are fully explored. As demonstrated in [|T^, ^ |TOl the LcWc term 
is particularly sensitive towards the r polarization, whereas the LeWe term 
determines the helicity of the r- neutrino. 

A detailed discussion of the matrix elements for the decay modes involv- 
ing different pseudoscalar mesons [/Cttttz/, KKKu, rjTfTiu] together with pre- 
dictions for the corresponding structure functions and angular distributions is 
presented in ||^. In this case, all 9 structure functions in (|^) are nonvanishing 
because of the interference of the anomaly with the axial vector contributions. 
The analysis of these distributions would allow to test the underlying hadronic 
physics, and to separate for example the contributions from the axial and the 
vector (Wess-Zumino anomaly) current. It is thus possible to confirm (not only 
qualitatively) the presence of the Wess-Zumino anomaly in the decay modes 
T -^ uKniT and r -^ uKittt. 

In IH^ the technique of the structure functions has been extended to the 
r -^ uuJTT decay mode. This allows to test the model for the hadronic matrix 
element, which involves both a vector and a second class axial vector current. 



2 Tau kinematics 

In all experimental and theoretical analysis of r decays it has been implicitly 
assumed that the r direction and its restframe cannot be fully reconstructed 




Figure 3: Kinematic configuration indicating the relative orientation of tfie 
liadronic tracks, tlie r directions and tlie vector dmm- 



and appropriately averaged distributions are considered. However, as shown in 
[p!7[] , impact parameter measurements allow to fully reconstruct the r direction. 
In events where both leptons decay semileptonically and where all hadron 
momenta are determined, the original r direction can be reconstructed up to 
a twofold ambiguity |1^, O, |15|, which can be resolved with the help of vertex 



detectors employed in present experiments. Several possibilities may arise: 

i.) If the beam spot is large compared to the typical impact parameter, 
the production vertex is unknown. Let us assume that both r decay into 
one charged hadron each and that both charged tracks can be measured with 
high precision. The direction dmin of the minimal distance between the two 
nonintersecting charged tracks (Fig.l) resolves the ambiguity and introduces 
two additional constraints that can be used to reduce the measurement errors. 
The r"*" and r" decay points and their original directions of flight are then 
uniquely determined. 

a.) Precise knowledge of the beam axis (corresponding to a beam spot 
of negligible size) leads to a further constraint resulting from the requirement 
that the reconstructed r axis and the beam axis intersect. If the production 
point would be known in addition then the momenta plus one charged track 
from the decay of only one r would allow to reconstruct the event and double 
semileptonic decays would lead to a large number of additional constraints. 



in.) If one (or both) r decay into several hadrons and if all momenta and 
two tracks (one from each side) are measured, the same reconstruction can be 
performed. 

Let us first consider case i) where all relevant aspects can be explained 
most clearly. The angles 6^ between the r^ and the hadron h^ directions 
respectively as defined in the lab frame are given by the energies of h'^ and h~ 

COS*' = y--~(l + ri)/2-, 



sin^^ 



(3\h'^x_ — r 



(1 - r2 )2/4 - (x_ - (1 + r^ )/2)V/3^ 



X- = Eh~/Er r_ = nih-lnir (11) 

and similarly for cos 6';^ and sin^;^. The velocity /3, and the boost factor 7 
refer to the r in the lab frame. 

The original r~ direction must therefore lie on the cone of opening angle 
9^ around the direction of h~ and on the cone of opening angle 9^ around 
the reflected direction of h'^ . The extremal situation where 9^ or 9^ assume 
the values or vr, or where the two cones touch in one line, lead to a unique 
solution for the r direction. In general a twofold ambiguity arises as is obvious 
from this geometric argument. The cosine of the relative azimuthal angle if 
between the directions of h^ and h~ denoted by n+ and n_ can be calculated 
from the momenta and energies of h^ and h~ as follows: In the coordinate 
frame (see Fig.l) with the z axis pointing along the direction of r~ and with 
n_ in the xz plane and positive x component 

/ sin( 
^- -n =\ I (12) 

cos^:^ 



\P-\ " ^ 



V+ _ ^ 

n+ 



I sin 9^ cos if 

sin 6^:^ sin yj | (13) 

V^ - cos 9^ 



\P+\ 
and cos(/? can be determined from 

n_n+ = - cos 9^ cos 9^ + sin 9^ sin 9^ cos 99 (14) 



The well-known twofold ambiguity in yj is evident from this formula. 

Additional information can be drawn from the precise determination of 
tracks close to the production point. Three-prong decays allow to reconstruct 
the decay vertex and the ambiguity can be trivially resolved. 

However, also single-prong events may serve this purpose. Let us first 
consider decays into one charged hadron on each side. Their tracks and in 
particular the vector dmin of closest approach (Fig.l) can be measured with 
the help of microvertex detectors. The vector pointing from the r_ to the r+ 
decay vertex 

d = T+-T- = -l\ (15) 

is oriented by definition into the negative z direction (/ > 0). The vector dmin 
can on the one hand be measured, on the other hand calculated from d, n_|_ 
and n_: 

dmin = d + [{dn+ n^n^ — dnJ)n_ (16) 

+ {dn^ n^n_ — dn^)n^]/{\ — {n_n^Y) (17) 

The sign of the projection of dmin on «+ x n^ then determines the sign of ip 
and hence resolves the ambiguity. 

dmin{n+ X n_) = I sin 9^ sin 9_ sin ip (18) 

The length of the projection determines I and hence provides a measurement 
of the lifetimes of r+ plus r_. Exploiting the fact that dn_ = —I cos 9^ and 
dn+ = I cos 9!^ the direction of d can be geometrically constructed by inverting 
®: 

d/l = dmin/l - [(cos 6*^ f?+n_ + cos6'^)n_+ (19) 

(— cos 9^ n^n_ — cos9^)n^]/{l — {n_n^Y) 

In addition two constraint equations may be derived by comparing dmin/l as 
calculated from the h+ and /i_ tracks with the direction calculated from (|T6|) 
with the help of 9^, 9^ and /. These might be used to constrain the events even 
in cases where initial state radiation distorts the simple kinematics described 
above. 

10 



As stated before, the locations of both r+ and r_ decay vertices in space 
are then fixed. If the beam axis is known with high precision (high compared 
to the decay length /) the lines between the two decay vertices and the beam 
axis intersect, providing one additional constraint. 

The generalization of this method to decays into multihadron states with 
one or several neutrals is straightforward: 6^^, 6^ and cos(p are fixed by the 
hadron momenta as stated above. Only one of the two solutions for ip is then 
compatible with dmin measured directly with the help of vertex detectors. 



3 T -polarization 

As discussed in chapter 1, the technique of structure functions serves as a 
usefuU tool to analyse the hadronic current and the r-polarization. Once the 
r restframe is reconstructed one may also proceed in a more direct way which 
allows to exploit maximally the spin information on an event by event basis. 
The current J^ depends generically on all hadron momenta. In [|T^ it has 
been demonstrated that the direction of the r-spin in the r-restframe in each 
event is given by 

h = H/u (20) 



with 



LU 



P^iU^-Ul); H = mr{fi''-fi) (21) 



n^ = 2[ir-N)j^ + iJ-N)j;-{r-j)N^]; 

Ul = 2\mt;''^j:j,N„. (22) 

The momenta P and N refer to the r and its neutrino. The V — A structure 
of the lepton current and a massless neutrino have been assumed for simplicity. 

For the decay into a single pion J^ oc Q^ and the direction of the vector h is 
identical to that of the momentum of the pion, its normalization equal to one. 
For the decay into several pions the current J^ and consequently the direction 
of h will depend on all pion momenta. It is, however, a simple exercise to 
demonstrate that 

P = 1 (23) 

also in the general case. This follows directly from 

(n^ - nj) (n'^ - n^^) = o (24) 

11 



and is a direct consequence of the V — A coupling at the lepton vertex and 
of the fact that no spin summation has been performed for the hadronic sys- 
tem. Also decay modes with kaons or rj^s only have maximal analysing power, 
corresponding to \h\ = 1. 

The argument is in particular also applicable for r -^ z/7ra;(— > Svr) if all 
four pion momenta are measured and the full matrix element including also 
the uj ^ Sn decay is used for the hadronic current. The argument would fail, 
if the spin of the u meson would be averaged. Correspondingly, |/i| 7^ 1 for 
leptonic decay modes. 

However, one important assumption has been made tacidly throughout: 
The form of the current must be known for the analysis of multi-hadron states. 
For 2 pions this form is essentially fixed, for 3 pions it is highly plausible. For 
multi-meson states, in particular those including kaons, a carefull test of the 
functional form of J is madatory. 

To summarize: Semileptonic decay modes of the r-lepton carry important 
information on hadron physics and r-properties. Measurements of the struc- 
ture functions are a convenient tool to constrain models and to determine 
the hadronic current in a model-independent way. Given sufficiently precise 
vertex detectors, double-semileptonic events can be fully reconstructed. Once 
the model for the hadronic current is reliable, the r spin can be measured in 
semileptonic decays on an event by event basis. 
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